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ABSTRACT 

Context. Although discovered 40 years ago, the emission mechanism responsible for the observed 
pulsar radiation remains unclear. However, the high-energy pulsed emission is usually explained 
in the framework of either the polar cap or the outer gap model. Here we explore an alternative 
model based on the striped wind. 

Aims. The purpose of this work is to study the pulsed component, that is the light-curves as 
well as the spectra of the high-energy emission, above 10 MeV, emanating from the striped wind 
model. Gamma rays are produced by scattering off the soft cosmic microwave background pho- 
tons on the ultrarelativistic leptons flowing in the current sheets. 

Methods. We compute the time-dependent inverse Compton emissivity of the wind, in the 
Thomson regime, by performing three-dimensional numerical integration in space over the whole 
striped wind. The phase-dependent spectral variability is then calculated as well as the change in 
pulse shape when going from the lowest to the highest energies. 

Results. Several light curves and spectra of inverse Compton radiation with phase resolved de- 
pendence are presented. We apply our model to the well-known gamma-ray pulsar Geminga. We 
are able to fit the EGRET spectra between 10 MeV and 10 GeV as well as the light curve above 
100 MeV with good accuracy. 

Conclusions. In the striped wind model, the pulses are a direct consequence of the relativistic 
beaming effect. It is a simple geometrical model capable to explain the very high-energy vari- 
ability of the phase-resolved spectrum as well as the shape of the associated pulses. Future con- 
frontations with observations at the highest energies and possibly polarization measurement will 
be important to discriminate between existing models. 

Key words. Pulsars: individual: Geminga - Radiation mechanisms: non-thermal - Gamma rays: 
observations - Gamma rays: theory - Stars: winds, outflows 

1. INTRODUCTION 



To date, seven pulsars are known to radi ate most of their energy in the gamma-ray band (about 
a few GeV), with high-confidence level dThompsorJ, l2003h . This high energy radiation is proba- 
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bly related to some ultrarelativistic particles which are flowing in the vicinity of the neutron star. 
Nevertheless, the detailed origin of the pulsed emission of a pulsar remains still poorly understood. 
A self-consistent description of the particle acceleration in the pulsar magnetosphere has also not 
been proposed yet. Many attempts have been suggested to explain the radiation process from a 
strong ly magnetized r otating neutron star . Among these, the polar caps dRuderman & Sutherland , 



19751) and outer gaps jCheng et aUll98q) are the most extensively studied models. The extremely 



stable time of arrival of the pulses supports the idea that the emission is probably coming from 
regions close to the neutron star surface. However, how close it should be is not clear yet. 

In the present w ork, we focus on an alternative model, the so-called striped wind model orig- 
inally introduced by ICoronitil (1 19901) and iMichell (119941) . In this picture, emission is expected to 



occur outside the light cyhnder but stifl in regions where the striped structure is well defined and 
rotates at the angular speed of the compact object. We use an explicit asymptotic solution for the 
large-scale field structure related to the oblique spUt monopole and valid for the case of an ultrarel- 
ativistic plasma (lBogovalovtll999l) . Due to r elativistic beamin g eff'ects, it has already been shown 
that this model gives rise to pulsed emission dKirk et al. , 2002 ) a nd can satisfactorily fit the optical 



polarization data from for instance the Crab pulsar (IPetri & Kirk , 



2005). 



It is usually claimed that the precise shape of the cut off" of the pulsed spectrum at the high- 
est energies (Ey > 10 MeV) will help to discri minate between diff'erent existing models. The 
outer gap scenario predicts an expo nential cut-off", Romanil ( 19961) whereas the polar cap model 



predicts a super exponential cut-off", iDaughertv & Hardina (119961) . In the near future, the phase- 
resolved spectra measured by GLAST will be a valua ble tool to confron t emission models to ob- 
servations at the highest energies (20 MeV-300 GeV) dThompsoni l2007h . The recent discovery of 
pulsed emission above 25 GeV in the Cra b will put some severe constraint on emission models 



dThe MAGIC CoUaboration: E. AUu , 



2008) 



Emission mechanisms outside the light-c y linder have a long history, as old as the discovery 
of the first pulsar Indeed, iMichel & Tuckeii (119691) proposed a supersonic magnetically domi- 
nated wi nd launched fro m the rotating neutron star, that eventually enters a shock and becomes 
subsonic. iLerchd ( Il970at) assumed that the low-frequency, large amplitude electromagnetic wave 



propagates in vacuum and is separated from the interstellar medium by an infinitely thin cur- 
rent sheet. This current sheet oscillates due to some magnetic pressure modu lation imposed by 

Lerchel Jl970biy) . 



the rotation of the pulsar. He also estimated the produced pulsed emission. In 
he extended this ide a to include finit e Larmor radius eff'ects which means a finite thickness of 
these current sheets. Dokuchaev et al 



( 119761) studied a similar phenomenon, namely the interac- 
tion between the spiral magnetic field structure and the interstellar medium in order to produce 
coherent emission. Coherent pulsed radiation emanating from a relativistically expan ding sheet 



Tademaru 



of cha r ges flowing along dipolar magn e tic field lines has been de eply investigated by 
( 1197 lb : iGrewing & HeintzmannI ( 1197 ih . iFerrari & Trussonil d 1975b found that pulsed radiation is 
possible at the light-cylinder and computed the associated synchro-c ompton s pectra . The forma- 
tion of neutral sheets and radiation in the wave zone is exposed in iMichell dl97lb . The pulsed 
high-energy emission could also be explained by magnetic reconnection (le ading to rel a tivisti c 
heating of particles) in the cur rent sheets just beyond the light-cylinder, see iLvubarskiil dl996h . 



Bogovalov & AharonianI d2000l) used inverse Compton scattering of the unshocked wind on low- 



frequency photons coming from the pulsar to compute the pulsed and unpulsed component of the 
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gamma ray emission. However, for some criti cism about pu lsed radiation coming from outside the 
light-cylinder and related beaming effects, see 



AronsI (119791) . 



Attempts to explain the pulse profiles and phase -resolved spectra of the Geminga pulsar have 
already been investigated bv lZhang & Cheng ( 2001 ) in the fra mework of a thick outer gap model. 
They used the data analysis performed bv lFierro et al.l (1 1998b according to EGRET observations. 
However, in their study, due to the geometry of the radiating region, no emission is expected in 
some phase of the rotation period of the neutron star. For instance, there is no leading wing before 
the first pulse or no trailing wing after the second pulse. They do not propose a fit of the light-curve 
at high energy. We will show that our model can explain both, the light-curve and the phase- 
resolved spectra. 

The paper is organized as follows. In Sec. |2] we d escribe the model , i.e. the magnetic field 
structure adopted, based on the asymptotic solution of iBogovalovl (119991) . The properties of the 
emitting particles in the striped wind, and the inverse Compton emission process are also discussed 
in this section. Next, in Sec.[3j we present the phase-resolved spectral variability and the change in 
the pulse profile of the light-curves with respect to photon energy. These results are then applied to 
fit the EGRET data {Ey > 10 MeV) for one gamma-ray pulsar, namely Geminga which is presum- 
ably not inside a synchrotron nebula. The conclusions and the possible extension are presented in 
Sec.g] 



2. THE STRIPED WIND MODEL 

The model used to compute the high-energy pulse shape and the phase-resolved spectrum aris- 
ing from the striped wind is briefly presented in this section. The geometrical configuration is 
as follows. The magnetized neutron star is rotating at an angular speed directed along the 
((9z)-axis i.e. the rotation axis is fl, - Q, e,. We use a Cartesian coordinate system with coor- 
dinates {x,y,z) and orthonormal basis The stellar magnetic moment m, assumed to be 
dipolar, makes an angle x with respect to the rotation axis, x - i^*7 m). This angle is therefore 
defined by cos;^' = m ■ e,/m where m - ||im|| is the magnitude of the stellar magnetic moment. 
Moreover, it rotates around c- at a speed f2* such that it is expressed as 

m - m [sin;if (cos(Q, f)ex + sin(Q, O^v) + cos;^'^,] (1) 

The inclination of the line of sight with respect to the rotational axis, and defined by the unit vector 
M, is denoted by ^, it Ues in the {Oyz) plane thus 

n - sin + cos^Cr (2) 

Moreover, the wind is expanding radially outwards at a velocity V close to the speed of light 
denoted by c. 

Our model involves some geometrical properties related to the magnetic field structure and 
some dynamical properties related to the emitting particles. Furthermore, in order to compute the 
light curves and the coiTesponding spectra, we need to know the emissivity of the wind due to 
inverse Compton scattering. This is explained in the next paragraphs. 
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2. 1 . Magnetic field structure 



The striped wind has been introduced bv lCoronitil (1 19901) and revisited by 



Michell (I1994 . Here, we 



ad opt a geometrical structure of the wind based on the asymptotic magnetic field solution given 
by iBogovalovl ( 119991) . This magnetic field model, although being only a crude approximation of 
the distant magnetic field, is nevertheless very convenient because it furnishes a simple analytical 
expression containing all the important features for an oblique rotator (relativistically expanding 
current sheets). The observed geometry is a well-known consequence of retardation effe cts winding 
up the stripes into a spiral. We emphasize that the striped wind solution as described in lBogovalov 
(119991) is not a necessary condition. A realistic model should lift some of the assumptions made, 
like a constant expansion speed for instance or infinitely thin current sheet. Nevertheless, the main 
characteristics of the pulsed emission properties would be preserved. However, in this first attempt 
to model the high-energy emission, we work with this simplified model. We leave improvements for 
future work. Thus, outside the light cylinder, the magnetic structure is replaced by two magnetic 
monopoles with equal and opposite intensity. The current sheet sustaining the magnetic polarity 
reversal arising in this solution, expressed in spherical coordinates (r, 6, (f) is defined by 



ct 

± arccos(- cot 6 cotx) H ip + 2ln 



(3) 



where (i - V/c, ri - c/Q.t is the radius of the light cylinder, t is the time as measured by a distant 
observer at rest, and I an integer Because of the ideal MHD assumption, this surface is frozen into 
the plasma and therefore moves also radially outwards at a speed V. Strictly speaking, the current 
sheets are infinitely thi n. However, as was a lready done for the study of the synchrotron polarization 
of the pulsed emission. iPetri & KirkI (l2005h . we release this restrictive and unphysical prescription. 
Indeed, the current sheet are assumed to have a given thickness, parameterized by the quantity A^. 
Moreover, inside the sheets, the particle number density is very high while the magnetic field is 
weak. In whole space, the magnetic field is purely toroidal and given by 

Rl 



(4) 



The strength of th e magnetic field at the light-cylinder is denoted by Bl. In the original work of 



Bogovalovl (1 1999b . the function 77^ is related to the Heaviside unit step function and can only have 



two values +1, leading to the discontinuity in magnetic field when crossing the surface defined 
by Eq. (|3]l. In order to make such transition more smooth, we redefine this unit step function by 
something going smoothly from -1 to +\ (and conversely) prescribed by the following expression 



77^ = tanh(A^(/r) 

i/r = COS co?,x + sin smx cos {(f - Q, (f - r/V)] 



(5) 
(6) 



With this formula, the transition layer has a thickness of approximately A^. The azimuthal depen- 
dence in (/; - Q, (f - r/V) represents am - 1 spiral wave propagating radially outwards at a speed 
V and rotating at the speed of the neutron star Q» . The structure of the azimuthal magnetic field in 
the equatorial plane is shown in Fig.[T] 



2.2. Particle distribution function 



The innermost regions of the pulsar magnetosphere is believed to be a site of high-energy pair 
production feeding the wind with ultrarelativistic electrons and positrons. For these emitting parti- 
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Fig. 1. Structure of the azimuthal magnetic field in the equatorial plane at a given time. Note the 
smooth polarity reversal when crossing the current sheets. The scaling by a factor r removes the 
radial decrease of the amplitude and helps to better visualize the thickness of the current sheets. 
Distances are normalized to the radius of the light cyhnder. 



cles, which are very dense in the cuiTent sheets where the magnetic pressure is weak, we adopt an 
isotropic distribution function in momentum space in the comoving frame of the wind. It is given 
by a power law in energy of index p, with a sharp low and high-energy cut-off, y,,,,-,, and jma^ re- 
spectively, such that the particle number density at time f and position r with energy between y and 

J + djis 

riAy, r, t) dy = Keir, t) j'P dy with 7,„,„ <y< y^ax (7) 

Ki.{r, t) is related to the number density of emitting particles in the cuiTent sheet. The particular 
form of the magnetic field in the current sheets, decreasing like a tangent hyperbolic function tanh 
as given in Eq. (|4]i, suggests to use a plasma density profile dictated by the exact solution of the 
relativistic Harris cu rrent sheet, namely a secant cosines hyperbolic function sech, see for instance 



Petri & KirkI (120071) . We thus adopt the following expression for the density 



^ (A^-A^o)sech2(A^<A)+A^o 
Ke(r, t) = (8) 

Nq sets the minimum particle density in the magnetized and cold plasma, i.e. between the current 
sheets, whereas defines the highest density inside the sheets. However, in order to allow different 
peak intensity in the light curves, we choose different maximum densities in two consecutive sheets, 
A'l and N2- This discrepancy could be explained by a different pair creation rate in the two polar 
caps, feeding the stripes with different numbers of particles. The radial motion of the wind at a fixed 
speed imposes an overall dependence on this quantity, due to the conservation of the number 
of particles. An example of particle density number, multiplied by this factor for the Geminga 
pulsar is shown in Fig.|2l Note, however, that adiabatic losses in the current sheets due to pressure 
work will cool down this distribution function in such a way that decreases like 1 / r^(/'+2)/3 ^}^^ 



19941) . As already done in a previous work, we assume the emission commences when the wind 



crosses the surface defined by r = ro > tl. Why should emission suddenly start at this arbitrary 
radius ro? If the stripes are responsible for the observed pulsed emission, it corresponds to regions 
where the wind is well form, i.e. outside the light-cylinder. Moreover, to explain the highly temporal 
stability of these pulses, emission should occur closest to the light cylinder because magnetic field 
is strongest there. Indeed, we do not expect strong radiation from far away regions because the 
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Fig. 2. The radial variation of the particle density number in the cold magnetized part (low density) 
and in the current sheets (high density) when looking in the equatorial plane for a given time. The 
scaling with the factor shows a constant density profile and a slightly different peak intensity in 
the current sheets. Distances are normahzed to the radius of the light cylinder. 



particle density number in the current sheets decreases like the square of the radius. Therefore, 
significant radiation is emanating from the base of the wind close to the launching region, probably 
near the light-cylinder Unfortunately, the transition region between the closed magnetosphere and 
the wind regime is poorly described, we only know that it should happen in the vicinity of the light- 
cylinder. The transition region would smoothly switch from an unpulsed component (in the closed 
magnetosphere) to a pulsed component (in the wind). In our model, we suppress this transition, 
making it of size zero. Thus, the relevant pulsed emission starts sharply when crossing the shell 
r - . Adding a smooth transition will not change the conclusions discussed in the next section. 
Note also that the whole wind is radiating, any plasma element located at a distance larger than 
ro will contribute to the total emissivity and not only when particles cross the shell tq. However, 
because of the fast decrease in intensity, only the first few stripes provide almost all the pulsed 
component. 

Let us now discuss the exact form of the emissivity functions in the striped wind for the inverse 
Compton radiation. 



2.3. Inverse Compton emissivity 

Let us first recall the general expression for Inverse Compton radiation of an electron embedded in 
an isotropic photon field. We assume an isotropic distribution of mono-energetic target photons s 
with density ny(e) in the observer frame. The number of photons of energy si scattered per unit 
time and per unit frequency by an ultra-relativistic electron is 

dN 3 o-jc 

^ = 57i/^^)«r©'^^ (9) 

For convenience, energies of photons are scaled to the rest mass energy of an electron c^ such 
th at £ = enipC^ and gi - S| c^. ctj is the Thomson cross section and the function / is defined 



m 



Blumenthal & GouldldlQTO) bv 



, (Ays of 

f{q) ^ l+q + 2q\nq-2q^+ {\ - q) (10) 

2(1 +Aysq) 

q = , (11) 
4ye(y- Si) 
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The power spectrum density of inverse Compton emission is therefore 

= (12) 
dtasi dtasi 

Summing the contribution from each lepton in the distribution function Eq. (|7]i, the total emissivity 
in the observer frame wiU be 

In the case of ultra-relativistic electrons, their rest mass energy is negUgible compared to their 
kinetic energy and thus they can be treated as massless particles. Their energy (or equivalently 
their Lorentz factor) then tran sforms according to the Doppler shift formula for photons y - 'Dj' 



dGeorganopoulos et all uOOll) where the Doppler factor associated with the relativistic expansion 



of the wind ^ = V/c with Lorentz factor F = 1 / ^1 - V^/c^ is given by 

£)= i (14) 

r(l-;8«) 

« being an unit vector along the line of sight joining the emitting element to the observer. Moreover, 
the associated distribution function expressed in the lab frame transforms according to 

nf\y)^D^n,(y/D)^K,ir,t)D"^^y-'' with Dy„in<y<Dy^, (15) 

For the remaining of this paper, we will only consider the Thomson regime for which ys<^ I. 
The last term on the RHS of Eq.lfTOli proportional to (1 - q) can thus be dropped. In this limit, the 
integration of Eq. (fTST l can be done analytically. To do this we need the following integral formula 

^'ic(-^) = J X^P'^ '''^ f{x)dx 

1 2x^ 2x Inx p-1 



(16) 



p + 1 p + 5 p + 3 {p + 3y 
The emissivity can therefore be written as 

jt'iei) = 2"-^ 3 cncnyis) D"^^ (-^J [F„(min(xi, 1)) - Fie(min(x2, 1))] (17) 

where the upper and lower limits are respectively xi - si/4-y^E and X2 - si /4 y^ s. the minimum 
and maximum Lorentz factor of the leptons are yj - Dy^i^ and 72 = i^Jmax- For frequencies well 
above the lower cut-off and well below the upper cut-off, the total intensity of inverse Compton 
scattering is then proportional to 

jt\sur, t) (X KAr, t) s^'"-'^^^ D^^^ riyis) (18) 



Comparing the power index of the Doppler factor in Eq. (3) of iPetri & KirkI (120051) and Eq. JTSl ), 
the dependence on p is different Ssync - (p + 3)/2 for synchrotron and Sic - (p + 2) for inverse 
Compton. For p > -I, beaming effect is weaker for synchrotron compared to inverse Compton 
because isync < ■Sic- Therefore, the pulses in the latter case will be more sharp than the pulses from 
the former. 
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2.4. Light curves and spectra 

Knowing the inverse Compton emissivity, the light curves are obtained by integration over the 
whole wind region. This wind is assumed to extend from a radius ro to an outer radius ro^t which 
can be interpreted as the location of the termination shock. However, because of the fast decrease of 
particle density with radius, there is no need to go such far in practice when performing the numer- 
ical integration. Moreover, the radius ro is interpreted as the transition region between the closed 
magnetosphere and the launch of the wind. We adopt a sharp unreaUstic transition by switching on 
the emission when crossing this surface. However, in reality we would expect emission also below 
this radius. But this region being poorly understood, we start with this crude model. The shape of 
the pulses will not be affected by this approximation, it is equivalent to neglect some contribution 
from the DC component. Therefore, the inverse Compton radiation at a fixed observer time t is 
given by 

I°^'(a),t) ^ j°^'(oj,r,t,^,)r^ sine drded(p (19) 

Jro Jo Jo 

The retarded time is expressed as fret = f - \\Ro - r\\/c ^ t - Rq/c + n ■ rjc. The approximation is 
vaUd if the observer, located at Rq, is very far away from the radiating system, » rout- Eq.(fT9]l 
is integrated numerically. We compute the inverse Compton intensity for several frequencies from 
far below the low cut-off frequency to far over the high frequency cut-off. We are therefore able to 
predict the phase resolved spectral variabiUty and the pulse shape simultaneously. The results and 
appUcations to the Geminga pulsar which is known to be a y-ray pulsars are discussed in the next 
section. 



3. APPLICATION TO y-RAY PULSARS 

We apply the aforementioned model to inverse Compton scattering of low-energy photons from 
the cosmic microwave background with typical energy of eqmb - Tqmb = 2.36 x 10 eV and 
energy density of 2.65 x 10^ eV/m^ by the ultrarelativistic leptons flowing in the current sheets of 
the striped wind. The CMB radiation field is assumed to be isotropic. We focus on the Geminga 
pulsar for which the useful parameters are listed in Table [U 

An important parameter of our model is the maximal Lorentz factor reached by pairs. Indeed, 
it directly determines the cut-off energy in the pulsed spectra. Note also that the minimal Lorentz 
factor does not play a significant role to fit the observations. We can find an order of magnitude 
of yinax in the following manner. Assuming that the highest energy emanates from the inverse 
Compton scattering with the CMB, we can deduce the maximal Lorentz factor of the electrons i n 



the observer frame. Indeed, no pulsed emission is seen above roughly 30 GeV ([Thompson , 



2004). 



The cut-off energy can be fixed at about 3 GeV. The maximum electron Lorentz factor is 
therefore obtained by 

r™t = (20) 

V 4 ^CMB 

Thus a typical value is y^^^ ^ 10^'^. The Lorentz factor of the wind, assumed to be constant. 



according to previous results iPetri & KirkI (120051) . can be set to a few tens. Adopting this value, 
the maximum Lorentz factor in the particle distribution function Eq. (|7]l in the comoving frame is 
typically ymax * 10^ 
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Pulsar 


P(ms) 


d (kpc) 


log(B) (T) 


Bi^ (T) (km) ( (degree) 


X (degree) 


Geminga 


237.09 


0.15 


8.21 


0.112 11312.4 90 


60 



Table 1. Properties of Geminga pulsar assu ming Rt - 10 km. Period, distance and magnetic field 



strength are taken from 



Taylor et al. 



(Il993l) . 



In our best fit, we choose an inclination of the magnetic moment with respect to the rotation 
axis of X - 60°. In order to obtain a phase separation of 0.5 between the two pulses, we have to 
adopt an inclination of the line of sight ^ - 90°. The Lorentz factor of the wind is F = 10. 

Results for the light-curve above 100 MeV and the definition of the different phase intervals 
are shown in Fig. [3] The rising and falling shape of both pulses are well fitted by our model. In one 
wavelength of the striped wind, or equivalently in one period of the pulsar, there are two current 
sheets which give rise to strong emission because of the high-density in it. This explains the double 
pulse structure which is observed at any frequency. Moreover, the two pulses are very similar in 
shape, with only a difference in maximum intensity. We explain this discrepancy by having slightly 
different particle density numbers within the current sheets. The corresponding phase-resolved 
spectra are shown in Fig.|4] The spectral variability is reproduced with satisfactory accuracy except 
for the OP phase for which the intensity is overestimated. 

The intensity and the width of both peaks is matched with good accuracy. Recall that the width 
is strongly influenced by the size of the current sheets and the Lorentz factor of the wind, due to 
relativistic beaming effects. For PI, the fitted spectrum lies within the error bars of the EGRET 
data whereas for P2, the slope is slightly overestimated. Nevertheless, the shape of the cut-off in 
both cases is retrieved. An excellent fit is also found for the trailing wing 1 (TWl) and the leading 
wing 2 (LW2). The bridge (BD) and the trailing wing 2 (TW2) follow the same trend. The cut-off 
although apparent in some phase intervals does not clearly show up in the off-pulse (OP) or in 
the leading wing 1 (LWl) interval. Observations form the future GLAST satellite will give useful 
information in this energy band and help to discriminate between different pulsar emission models. 



. 8 
I 0.6 



I P2 I TW2 




I LW2 



Fig. 3. Gamma-ray light curve above 100 MeV of Geminga fitted with the inverse Compton emis- 
sion from the striped wind. 



4. CONCLUSION 

In the striped wind model, the pulsed high-energy emission from pulsars arises from regions well 
outside the light-cylinder The geometry of the wind does not depend on the precise magnetospheric 
structure close to the neutron star surface. Therefore, a detailed model of the physics in the inner 
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Fig. 4. Phase-resolved inverse Compton emission from the Geminga pulsar for different phase inter- 
vals: bridge (BD), leading wing 1/2 (LW1/LW2), off -pulse (OP), peak 1/ 2 (P1/P2), trailing w ing 1/2 
(TW1/TW2). Exact definition of these intervals and data are taken from lpierro et al. (1998). 



pulsar magnetosphere is not required to make robust predictions concerning radiation emanating 
from the striped wind. By computing the inverse Compton emission on the CMB photons in the 
Thomson regime, we were able to fit the EGRET data of the light-curves as well as the phase- 
resolved spectra for at least one gamma-ray pulsar, like Geminga. 

Measuring the polarization properties of this high-energy emission would provide a valuable 
diagnostic of the physical process responsible for these MeV-GeV photons. In the future, experi- 
ments are planned to measure the polarization in X-rays or gamma-rays, like the PoGOLite satellite 



jKamae et al. 



,|2007D. 

Because the radiation is mostly emanating from the basis of the striped wind where the density 
of particle is highest, a better description of the transition between the closed magnetosphere and 
the wind is required to obtain a detailed view of the magnetic field configuration. 

Note also that the same kind of computation can be performed for the Crab pulsar for which 
we know the polarization properties in optical bands and can therefore predict the polarization at 
higher energies by extrapolation. 
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